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The effect of interference stabilization is shown to exist in a system of two atomic levels coupled 
by a strong two-color laser field, the two frequencies of which are close to a two-photon Raman- 
type resonance between the chosen levels, with open channels of one-photon ionization from both of 
them. We suggest an experiment, in which a rather significant (up to 90 %) suppression of ionization 
can take place and which demonstrates explicitly the interference origin of stabilization. Specific 
calculations are made for H and He atoms and optimal parameters of a two-color field are found. 
The physics of the effect and its relation with such well-known phenomena as LICS and population 
trapping in a three-level system are discussed. 

PACS numbers: 



I. INTRODUCTION 

A. Interference stabilization 

Interference stabilization of Rydberg atoms, or strong- 
field suppression of photoionization, is known I'], '3| to 
be a phenomenon related to the coherent re-population of 
levels neighboring to the initially populated one. Such a 
re-population arises owing to Raman-type transitions via 
the continuum and in the case of a single-color field it can 
be efficient only if the field is strong enough. Specifically, 
the strong-field criterion for the effect of interference sta- 
bilization is formulated qualitatively as the requirement 

(n) 

for the ionization width of the initially populated 
atomic level En to be larger than the spacing between 
neighboring levels. 



(1) 



where n is the principal quantum number. The ioniza- 
tion width is determined here as the rate of ionization 
calculated with the help of the Fermi Golden Rule. How- 
ever, repopulation of neighboring Rydberg levels is pro- 
vided, actually, by the off-diagonal terms of the tensor of 
ionization widths Fj-" '"^ In the approximation of adi- 
abatic elimination of the continuum (which includes, in 
particular, the well-know rotating wave approximation, 
[J ) this tensor is determined as a direct generalization of 



the Fermi Golden Rule expression for F, 



(») 
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where Sq and u are the laser field-strength amplitude and 
frequency, d is the projection of the atomic dipole mo- 
ment upon the direction of light polarization, E is the en- 
ergy of the atomic electron in the continuum, and atomic 
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units are used throughout the paper if not indicated dif- 
ferently. So, the next crucial assumption in the theory 
of interference stabilization is that all the components of 
the tensor Q are approximately equal to each other 



F, 



(n' ,n) 



(3) 



This assumption is pretty well fulfilled for high atomic 
Rydberg levels, n, ri' 3> 1, \n — n'\ n (see explanations 
in 0). It should be noted also that for Rydberg levels 
their ac Stark shift, as well as the shift of the ionization 
threshold are equal approximately to the ponderomotive 
energy el/4:Uj'^ and identical to each other, and this com- 
mon shift does not affect either the dynamics of photoion- 
ization from Rydberg levels or the effect of interference 
stabilization. 

The simplest model, in which the effect of interference 
stabilization exists, is the model of two close atomic levels 
El and E2 connected with each other by the Raman-type 
transitions via the continuum, for which the conditions 
and (O are fulfilled and the ac Stark shift has the 
same features as described above for Rydberg levels and, 
actually, can be ignored. 

Both in two-level and multilevel systems there are sev- 
eral different theoretical approaches one can use to solve 
the problems of strong-field photoionization and stabi- 
lization. One of them is based on the use of quasienergy 
or " dressed-state" analysis. The total wave function of 
an atomic electron in a light field can be expanded in a 
series of the field-free atomic eigenfunctions 



5* = C„(t) ipn + continuum. 



(4) 



In the approximation of adiabatic elimination of the con- 
tinuum equations for the coefficients C„ (t) are stationary, 
and in the simplest case of the two- level system they have 
the form 



tCi (t) - Si Ci (i) = - f F [C\ (t) + C2 it)] 



(5) 



mckel^ahETrur 



ic\ it) - E2Ciit) ^ - [Ci it) + C2 it)] , 

where the approximation Q is assumed to be fulfilled. 
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As equations ^ are stationary, they have solutions 
of the form Ci^2 k exp(— «7i), where 7 is a complex 
quasienergy. When this exponential dependence on t is 
substituted into Eqs. they turn into a set of two 

algebraic homogeneous equations, which has a nonzero 
solution if its determinant turns zero. This is the condi- 
tion from which the two quasienergies of the field-driven 
two-level system have to be found, and the result is given 
by 

7± = 5 {^1 + £^2 - i r ± ^iE2 -Ei)^-r^Y (6) 

From here we see that, indeed, a drastic change in the 
form of the solutions occurs when the interaction con- 
stant r becomes larger than the level spacing E2 — Ei. 
The point T — E2 — Ei is the branching point, below 
which (at T < E2 — E1) the root square is real in Eq. ||HJ), 
whereas above the branching point (at F > £'2— -Ei) it be- 
comes imaginary. The imaginary parts of the quasiener- 
gies 7± © are shown in Fig. 1, and they determine 
the field-dependent widths of the two quasienergy levels 
r±(eo) = 2|Im[7±(F)]|, where F oc ©. One of the two 
branches arising dX T > E2 — Ei (7+(F)) corresponds 
to a narrowing quasienergy level whose width F_|_(£o) 
falls with a growing field-strength amplitude. This corre- 
sponds to an increasing life-time of this quasienergy level 
and to stabilization of an atomic population at this level. 



Im(yJ 




Fig. 1. The functions Im[7±(r)] ©. 



B. Laser-Induced Continuum Structures, 
Autoionizing Resonances, Dark States, and 
Population Trapping 

Though rather attractive by its simplicity, an isolated 
two-level system obeying the requirement © hardly can 
be easily found in the usual atomic spectra. This is the 
reason why here we consider another scheme, in which 
two atomic levels with significantly different energies and 
bound-free dipole matrix are connected with each other 
by Raman-type transitions via the continuum in a two- 
color field (Fig. 2). Such a scheme has been widely dis- 
cussed in literature 

connection with the phenomenon of Light-Induced Con- 
tinuum Structure (LICS), briefly outlined below. The 



process we suggest and investigate, as well as its similar- 
ity and differences with LICS are discussed in Subsection 
C. 

In LICS, one of the two fields of a two-color light is as- 
sumed to be strong [(w2, £2), the pump] and the other one 
- weak [{loi, ei), the probe], where wi,2 and £12 are the 
corresponding frequencies and field-strength amplitudes. 
In a scheme of Fig. 2 A is the Raman-type resonance 
detuning. 

/^^E2+UJ2-Ei-UJi (7) 
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Fig. 2. A scheme of two atomic levels under the conditions 
of a Raman-type resonance in a two-color filed. 

Under the action of the pump the first Floquet satellite 
of the level E2 takes the form of an autoionizing- like level 
E2 -\- L1J2 at the background of the continuum with the 

(21 

width equal to the ionization width of the level i52, r,- = 

Fp'^-* |5J with Eo substituted by £2- If all the population 
is concentrated initially at the level Ei, the probe field 
ionizes the atom and takes the Floquet satellite E2 + i02 
for an almost real autoionizing level, which gives rise to 
the typical asymmetric Fano-profile-like structure of the 
dispersion curve 'Wi{ijJi) (Fig. 3). 




E2+(i>2-Ei 



Fig. 3. The Fano profile of the dispersion curve 'Wi{u\) at an 
autoionizing-like resonance. 
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The Fano minimum of the curve Wi{uji) arises owing to 
interference of direct and indirect transitions to the con- 
tinuum {El E and Ei ^ E' ^ E2 ^ E). But 
this is not yet a stabihzation understood as an increasing 
suppression of ionization with a growing Hght intensity. 
An experiment we suggest and discuss below can demon- 
strate expHcitly such an interference suppression and sta- 
bihzation of an atom in its bound states. We wiU assume 
that both fields can be equally strong and in this sense 
the effect we consider can be referred to as a strong-field 
Lies. 

It should be noted that the channel of ionization shown 
in Fig. 2 by a dashed line (taken into account for the first 
time in Ref. 0) determines a nonzero height of the Fano 
curve in its minimum. As we will see below, this is just 
the competition of this non-interfering channel of ioniza- 
tion with interfering ones that is crucially important for 
optimization of stabilization in its dependence on light 
intensities. 

At last, another well-known analogue of LICS is re- 
lated to single- and double strong-field resonances at real 
autoionizing atomic levels 7, 14, 15, lilli^. The physics 
of this phenomenon and LICS are close though there 
are evident differences concerning mechanisms of level 
broadening. These differences make autoionizing reso- 
nances not as closely related to the phenomenon under 
consideration as LICS. The same can be said about dark 
states and population trapping in a three-level system 
[H m 113, ij. The physics of all these phenomena 
is alike though important details are different. In par- 
ticular, this concerns intensity-dependent mechanisms of 
level broadening and an important role of the noninter- 
fering ionization channel specific for the scheme under 
consideration and missing in a three-level scheme. Be- 
sides, the atomic continuous spectrum is so much wider 
than any discrete third level that this makes character- 
istic intensities in the phenomenon to be discussed ab- 
solutely different from those in the population trapping 
effect. 



C. An experiment we suggest 

Let us assume that the initially populated level in a 
scheme of Fig. 2 is At the first stage, by considering 
ionization of an atom by the field Si alone (with £2 = 0), 
we select the peak light intensity Ii = el/Sn and pulse 
duration r high and long enough to provide almost com- 
plete ionization of an atom by a pulse, Wi^i, /2 = 0) = 1. 
Then, by adding the field £2 we expect that under proper 
conditions, owing to interference effect, the combined ac- 
tion of two fields will result in a significant suppression 
of ionization. In other words, we expect that the func- 
tion Wi{Ii, I2) in its dependence on I2 at a given se- 
lected value of Ii will start from 1 at /2 = 0, then it 
will have a minimum at some intensity I2 ~ l2a ^ Ii, 
and then, at higher intensity I2, Wi{Ii, I2) will return to 
one again. The region around /20 will be interpreted as 



a stabilization window. The interference origin of stabi- 
lization is evident, because with a growing I2 (at a given 
/i), we increase an energy that can be put into an atom. 
But, counterintuitively, this additional energy results in 
a slower rather than faster ionization, and this can be 
explained only by interference effects. 

We will find values of the resonance detuning A and the 
ratio of intensities I2 / Ii which optimize the stabilization 
effect. 



II. THE MAIN EQUATIONS 

Compared to a simplified system of two close levels in 
a single-color field described in the Introduction, to char- 
acterize appropriately a system of Fig. 2 in a two-color 
arbitrary strong field, in addition to ionization broaden- 
ing and mixing of levels, we have to take into account 
also their shifts and mixing arising owing to the ac Stark 
effect. Both these effects can be described in terms of the 
complex polarizability tensor a^.j, i,j — 1,2, 



= fdE \d,E\' (t^^ ^ + -W-4r—) (8) 

J \E — Ei — uj — 10 E — Ei + uj J 

and 

a2i= [ dE d2EdEi[ -p, tt— 

J \E — El — uji — 10 

where integration over E includes summation over inter- 
mediate discrete states. 

For the two-color scheme (Fig. 2), similarly to Q, in 
the rotating wave approximation the wave function of an 
atomic electron can be written as 

^ = Ci{t) e^'^^Vi + C2{t) e'"^*?A2 + continuum. (10) 

As well as Eqs. (O, equations for the probability ampli- 
tudes Ci [t) and C2 (i) are obtained from the Schrodinger 
equation with the help of the procedure of adiabatic elim- 
ination of the continuum. In terms of a^j these 
equations can be presented in the form 

iCi - {Ei{t) + LUi)Ci = -\eio{t)e2Q{t)ai2 C2, 

(11) 

iC2 - {E2{t) + UJ2)C2 - -i£io(i)£2o(i)a2i Ci, 

where Eiit) are the slowly time-dependent adiabatic 
complex energies of the ac- Stark-shifted and broadened 
levels 

E,{t) =E,- \[a,{LOi)elS) + a,(cc;2)£loW}- (12) 
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III. QUASIENERGIES 



IV. PROBABILITY OF IONIZATION 



For the time-dependent pulse envelopes Sio{t) and 
£2o(0; Eqs. have to be solved as the initial- 

value problem. In the model of constant field-strength 
amplitudes these equations have stationary solutions 
Ci^2 oc exp(— «7t), where, as previously, 7 is the complex 
quasienergy for which we get the solutions generalizing 
those of Eq. ® 



where 



(13) 



(14) 



and A is the time-dependent complex detuning for the 
ac-Stark-shifted and broadened levels ifT^ 



A = E2 + LiJ2 — El ~ UJl. 



(15) 



Imaginary parts of the energies Ei (I12|l are related 
to the ionization widths determined by imaginary parts 
Im(ai) = a" of the polarizabilities ai (|SJ) 



Im(Si) = -iri, 

im{E2) = -\ (r^'^ + r. 



(2)^ 
2 ) ' 



(16) 



where 



Ti = ia'/(c^i)£?o, 



-(1) - 1^,"^ 



2 - 5"2('^i)£i0' and ^ = ia^'(u;2)eio• 



(17) 



^ 2 

The width r2^'' is determined by transitions from the level 



E2 under the action of the field £10 (the dashed line in 
Fig. 2). As E2 > El and lui > lu2, typically, 



a'^icji) < a'licji) or F^^' <C Fi. 



(18) 



Similarly to H17|l . the off-diagonal component of the po- 
larizability tensor ||^ determines the off-diagonal compo- 
nent of the ionization-width tensor 



F12 — ^Q;"2£io£2 — v riFj^^ 



(19) 



which assumes, in particular, that ai2 — 

Imaginary parts of quasienergies 7± (|13|l are related to 
the width of quasienergy levels T± 



-2Ini(7i 



(20) 



The described above quasienergy solutions of Eqs. Ijllll 
are most appropriate for solving the initial- value problem 
in the case of pulses with rectangular envelopes £1,2 o(Oi 
with sudden turn-on and turn-off (at t — and t — t) and 
£1,20(0 = const at < i < T. With known quasiener- 
gies 7± (|13|l the time-dependent probability amplitudes 
(^1,2(0 to find an atom in its bound states V'l and ip2 can 
be presented in the form 

Ci.2{t) = A+^expi-i-f+t) + Al^expi-i-f^t), (21) 

where Af2 are constants to be found from the initial 
conditions 



A^ 



a; 



1 and Ao 



(22) 



and from equations connecting A^ with The latter 
follow, e.g., from the first of Eqs. (|ll(l 

j±Af - {El + uji)Af = -i£io£2oai2 Af. (23) 

The total residual probability Wres (t) to find an atom 
in bound states at t = r is given by the sum of partial 
probabilities wi{t) and 1^2(7) to find an atom at levels 
El and E2, 



Wres(j) = Wi{t) + W2{t), 

where 

Wi,2{t) = |Ci,2(t)P 



(24) 



-I 7+ r 



A, 26 



-I 7_ r 



(25) 

The probability of ionization is given by Wiir) = 1 — 
Wresir). Eqs. and are easily solved to give, 

explicitly. 




and 



W2 (r) 



12^10^20 
16 £12 



(26) 



(27) 



In the case of pulses with smooth envelopes £1,2(0 
quasienergy solutions are not so useful for solving the 
initial-value problem, and one has to solve directly 
Eqs. (|lll) for the time-dependent probability amplitudes 

Cl,2it). 



V. SCALING EFFECT AND RELATIVE UNITS 

With the help of a phase transformation 

C,{t) = exp {-i{Ei + uj)t} Ai{t) (28) 
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equations can be reduced to an asymmetric form 



iAi 



and 



i{«iK)e?oW+ai(w2)eioW}^i 
= -jSio{t)e2o{t)ai2 Ai 

= -■|eio(i)e2o(0"2i ^i, 



(29) 

where, as previously, A is the weak-field two-photon- 
resonance detuning iQ. 

Though not as nice as Eqs. 1)291) are more con- 

venient to describe the scaling effect existing in the sys- 
tem under consideration. Let us assume that both pulse 
envelopes eio(i) and £2o(0 depend on time t only via 
the ratio t/r, where r is the pulse duration common for 
both high- and low- frequency pulses. Then, evidently, 
the arguments of the functions eio(^) and e2o(i) do not 
change if we divide both t and r by the same factor A, 
t ~^t/\ and T r/A. Moreover, one can see easily that 
Eqs. (|29|l do not change too if we multiply simultane- 
ously both low- and high-frequency pulse peak intensities 
^1,2 = cef 20/8'"' ^iid the weak-field detuning A by the 
same factor A. So, the solutions of Eqs. (|29|) are invariant 
with respect to the scaling transformation: 

A ^ AA, £2 2„ ^ A£?,2o, ^ ^ t ^ i/A. (30) 

with an arbitrary A. This scaling effect can be important 
for practice: parameters of an assumed experiment can 
be varied to choose the most convenient conditions for 
observation the two-color stabilization effect discussed 
in this paper. In particular, by making laser pulses 
longer, one can use rather moderate-intensity lasers, as 
it's shown below. 

Owing to the the described scaling effect, it is con- 
venient to introduce and use the dimensionless ratio of 
intensities x, interaction time 9 and detuning 6, 



h A 

X = — , e = T ■ h, S = —, 

dimensionless complex quasienergies 



y± 



1±- El- LUl 

h 



(31) 



(32) 



detuning for the ac Srark shifted and broadened levels 



= 5 - i|Q;2(ti^i) -I- ai{i^i) + [Q:i(ti'2) + 0^2(^2)] a;|, (33) 
and widths of the fully dressed quasienergy levels 



; [q" (tiJi) + a{{uj2)x\ - 5" =F Im 



1 2 ' 

I GL 1 n X 



(34) 

where Ii^ 2, A, r, 7-1-, and T± are in atomic units. Defined 
in such a way, quasienergies y± and widths g± depend 
only on two parameters, x and (5, whereas the probabil- 
ity of ionization Wi and the residual probability to find 
an atom in bound states Wres H24|) - (|27|) depend on three 
parameters, x, 5, and 9. 



VI. PULSE SHAPE 

The concepts of quasienergies and quasienergy func- 
tions are very fruitful for an analysis exploiting a model 
of a rectangular pulse envelope. Such an analysis is useful 
for clarification of physics of the phenomenon under con- 
sideration. However, more realistic laser pulse shapes are 
characterized by smooth envelopes. To investigate a sen- 
sitivity of the results to be derived on the pulse shape and 
its smoothing, we will consider pulse envelopes £1.2(0 of 
the form 



(1 + a) sin^ 



£1,2 a(t) = £1,2 X 



7r(iV(a)--^- 



l + a sin 



(35) 

where £1,20 = const., -l/2N{a) < t/r < l/2iV(a), 
N{a) is the normalization factor 



such that 



-/2N 



-t/2N 



'1,2 



2 + 3a 
2(l + a)3/2 



o(i) dt^TX £? 2 0: 



(36) 



(37) 



and a is a smoothing parameter. At a — > 00, N{a) — >1 
and the envelopes £1. 2 o{t) (|35|l turn into the rectangular 
ones. At a = 0, A^(0) = 3/8 and £1,2 o(i) ES) takes the 
form of a pure sin^ pulse envelope 



£1,2 o(i) = £1,2 X sin^ 



3t 1 



(38) 



g± 



-2Im[2/j 



By definition, for the pulse envelopes of the form 1)35)1 . 
at all values of the smoothing factor a the peak values 
of the field strengths and areas under the curves ef 2 i^) 
(field energy per unit cross section) are kept constant and 
equal to £1,20 and £120^'''' correspondingly (see Fig. 
4). Hence, at all a the parameter t can be interpreted 
as the pulse duration determined by the condition 137)1 . 
For all a the dimensionless pulse duration is determined 
as previously, 9 = Ii x t with both Ii and r taken in 
atomic units. 
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Fig. 4. Pulse envelopes 13511 at three different values of the 
smoothing factor a. 



VII. OPTIMIZATION OF THE 
LEVEL-NARROWING EFFECT 

Optimization of the stabilization and level-narrowing 
effects assumes minimization of the width of the narrower 
quasienergy level 5+ (5, x) with respect to two variables, b 
and X. As a first step in a solution of this problem, let us 
minimize 5+(<^, x) with respect to the field- free detuning 8 
at a given x. Typically, at any given x, in dependence on 
(5, the functions g±(5, x — const) have well pronounced 
minimum and maximum (Fig. 5). 




Fig. 5. The functions g±{S) = g±{5,x = cowat') II33II . 



Specifically, the curves of Fig. 5 are calculated for a 
He atom at x — 0.1. Details of these and many other 
calculations, as well as the data about frequencies, atomic 
levels, and polarizability tensors are given in the following 
Section. Here the picture of Fig. 5 is shown as a typical 
example of the dependencies g± {5, x — const) . To find 
a position of the extremum shown in Fig. 5, Soptix), 
we have to solve the equation dg±{5)/dd = 0. Direct 
calculations show that this condition is satisfied if Im{ [S ■ 
a*2]^} = 0, which gives two equations: Re[6 ■ 0*3] — 
and lm[6 ■ a^j] — 0- It can be checked directly that 
only the second of these two equations corresponds to 
the extremum we are looking for, and this equation gives 



^'12 S" - a'1'2 S' = 0, or S,pt = ^ <5 ", (39) 

"12 



or 



SoAx) = - o^'ii^i) - ^K'K) - «'i'(^i)]} 



+ j{a',{u;2) - a[iu;2) - ^^(0.2) - a'liu;^)]}. (40) 
* ^ "12 ■' 

As explained above, Sopt{x) is a value of the field-free de- 
tuning, at which the width of one of the two quasienergy 
levels (7+) has a minimum with respect to S at arbitrary 
given X. As it's seen from Eq. 140(1 the optimal detun- 
ing Sopt{x) is a linear function of the ratio of intensities 

X = l2/Il- 

The second step in optimization conditions for stabi- 
lization requires minimization of the width calculated 
at (5 = 6opt{x) with respect to the variable x. With 
the help of the second equation H39|l the "5-optimized" 
widths g±{Soptix),x) can be reduced to the following 
rather simple form 

9±{Sopt{x),x) = 



ia'/(^i)-<5;;,(x) 



[S'^x)]' 



1 // 2 



where 



1 / 



Soptix) = -^[a^K^i) a'i'(c.i) + a'2'(c.2)x 



(41) 



(42) 



and in Eqs. H41|) and H42|) we have put a"(w2) = 0, which 
is true in the case lu2 < l^'il- 

By using Eqs. and we can find easy the 

asymptotic expansion of the function g^{Soptix),x) in 
powers of l/x at large x, x 3> 1. The constant term 
in this expansion vanishes and the first non-zero term is 
given by 



g+{Sopt{x),x) 



1 ai"{LJi)a2"iuJi) 



2x 



012 "(t^2) 



(43) 



This result shows that at large x the width of the nar- 
rower quasienergy level 5+ decreases and tends zero as 
l/x. This is an indication of a possibility of the un- 
limited narrowing of this quasienergy level and, hence, 
achievement of a very high degree of stabilization. Real 
limitations of narrowing and stabilization are determined 
by the applicability conditions of the model. E.g., at very 
large values of x the second-field intensity can become too 
high for the ATI processes in this field to be ignored. In 
accordance with Eq. (|40|l . at the optimal conditions the 
increase of x is accompanied by a linear increase of the 
field-free detuning 5. This gives other limitations for the 
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growth of x: at sufficiently large detunings the influence 
of atomic levels different from Ei and E2 and not taken 
into account in the model can become important. At 
last at very large S even the used above rotating wave 
approximation can become invalid. But, on the other 
hand, these limitations are not too severe, and rather sig- 
nificant level of narrowing and stabilization of an atom 
can be reached under quite realistic conditions. These 
conclusions, as well as the general result about asymp- 
totic decrease of the narrower-level width at large x are 
confirmed and specified in the following Section by direct 
numerical calculations for Hydrogen and Helium atoms. 

It should be noted that the curve Sopt{x) includes the 
point {So, xo) where the complex detuning between the 
ac-Stark shifted and broadened levels (|f 5|l turns zero, 
(5 = 0. At this point 



Xq 



(44) 



0^2(^2) - «! (W2) 

and, in accordance with Eq. H40|l . 
^0 ^ ^{a'^iuJi) - a[{uJi) + xo[a'2{uj2) - a[{LJ2)]} ■ (45) 



VIII. NUMERICAL CALCULATIONS FOR A 
He ATOM 

All the required numerical data are known for selected 
levels in Hydrogen and Helium atoms. To a very large 
extent, the results obtained for these two atoms appear 
to be very similar. For this reason, below, the results of 
calculation of a Helium atom are described in details, and 
then a short sketch of calculations for a Hydrogen atom is 
given to demonstrate mainly the arising differences and 
peculiarities of each atom. 



The relative width of the fully dressed quasienergy lev- 
els g± H34|l are shown in Fig. 6 in their dependence on x 
at S = 6q. 




0.1 0.2 0.3 0.4 0.5 0.6 

Fig. 6. The relative widths of quasienergy levels of a He atom 
g±{x) calculated a.t S = So 14711 . 

At this value of the detuning 6, the curves g+{x) and 
g- (x) cross each other twice. The left and right crossings 
turn into the avoided crossings, correspondingly, sd S > 
— 11.3 and S < — 22.4, and both crossings never turn into 
avoided crossings together. 

For the polarizability tensor of Eq. H4t)|) , the expression 
(|40|l for the optimal- narrowing detuning 5opt{x) takes the 
form 



-0.26 - 83.572;. 



(48) 



In the picture of Fig. 7 the width g+{5, x) of the nar- 
rower quasienergy level is plotted in its dependence on 
the intensity ratio x at three different given values of the 
field- free detuning 5. The curve at 5 = —II differs 



A. Widths of quasienergy levels 

Let us consider the following two levels of a He atom 
and the following two frequencies: Is2s = Ei and Is4s = 
E2 and wi = 8.44 eV (the second harmonic of a dye-laser) 
and lj02 = 1,17 eV (Nd:YAG laser). For these levels and 
frequencies all the polarizability tensor components lO, 
are known [l^,!!^, and in atomic units they are equal 
to 

ai(wi) = -30.42 -hi 22.65, ^1(^2) = -236.6, 
a2{uJi) = - 45.66 -I- i 3.21, ^2(^2) = - 479.96 -I- i 124.55, 
and ai2 = 38.74 -f- i 53.07. 

(46) 

The point where the dressed-level detuning A (|I5f) (or 
S (|33|l ) turns zero, A (5 = 0, has the following coordi- 
nates in the plane {x, 6}: 

Xo = 0.156 and = - 13.3. (47) 




X 



0.1 0.2 0.3 0.4 



Fig. 7. Three typical curves g+{S = const, x) for He calcu- 
lated at 5 = —11, —15, and — 23 (from the left to the right) 

qualitatively from two other curves. The difference arises 
because for this curve the detuning S is large enough 
for the left crossing of Fig. 6 to turn into the avoided 
crossing. The sharp peak of the curve g+{5 — —11, x) at 



8 



Fig. 7 is an indication of the root-square branch-point- 
like behavior, which takes place at (5 = —11.3 and which 
is similar to the root-square branch-point behavior of the 
curve of Fig. 1 for 7^ in an idealized two-level system of 
Section lA. 

The minima of all three curves at Fig. 7 are seen 
to be getting the deeper the larger is \d\. Positions 
of these minima, x™™((5), are determined by the opti- 
mal detuning (|48(l and can be found from the equation 
S = Soptix). The X— dependent width of the narrower 
quasienergy level, minimized with respect to 6, is given 
by f/™*"(a:) = g+{Soptix), x), and its dependence on x is 
given by the curve of Fig. 8. 




_l I I I I 

0.5 1.0 1.5 2.0 2.5 3.0 



Fig. 8. The width of the narrower quasienergy level of a 
Helium atom minimized with respect to the detuning 

S, in its dependence on x. 



(|^ , , (j?7|) , and the results of the corresponding cal- 
culations are shown in Fig. 9. The three resonance-like 



w„,(8, x) 

1.0 - 




_I I I I L_ 

2 4 6 8 10 



Fig. 9. The residual probability to find a He atom in bound 
states calculated at 5 — —100, —250, and —500 (from the left 
to the right) and 6 = Sopt{x) 1481 1 (the upper curve), 9 = 0.1. 

curves correspond to three different values of the field-free 
detuning 5. The envelope of peaks of these curves is the 
maximized residual probability equal to Wresi^optix), x) 
with Sopt{x) given by Eq. (|^ . These results show 
that the function WresiSoptix),x) monotonously grows 
approaching one at large values of x. This means that 
under optimal conditions stabilization of a He atom in a 
two-color field can be very high, more than 90%. 

The picture of Fig. 10 shows the distribution of the 
residual probability between the levels Ei and E2 ■ 



A monotonous fall of the function g™"^ (x) means that 
the x— dependent width of the narrower quasienergy 
level, minimized with respect to the detuning S, tends 
to zero at asymptotically large values of x. In other 
words, in the framework of the used model one of the 
two quasienergy levels of the system can be narrowed un- 
limitedly by means of increasing the ratio of intensities 
X = I2/I1 and the field- free detuning \6\ in such a way 
that the equation d — 5 opt {x) remains satisfied. Real lim- 
itations of such a narrowing are determined only by the 
model applicability conditions: (i) at very large values 
of x the second-field intensity I2 will become too high to 
ignore above-threshold ionization produced by this field 
and (ii) at very large values of \5\ other levels ignored 
above can become more important and, at |A| ~ ^2, the 
rotating-wave approximation can become invalid. But at 
sufficiently long pulse durations r and low first-field in- 
tensity Ii these limitations are not too severe, and the 
achievable degree of narrowing can be rather high. 



B. Probabilities of ionization and non-ionization 

The residual probabilities to find a He atom in its 
bound states after interaction with a two-color field in 
the case of a rectangular envelope is determined by Eqs. 




2 4 6 8 



Fig. 10. The ratio ■W2/W1 for a He atom a,t S — 5opt{x) 1481 1. 
e = 0.1. 

Under the conditions of optimal stabilization (S = 
^optix)) at sufficiently high value of the intensity ratio 
X, the ratio of probabilities Wi/w2 falls tending asymp- 
totically to zero. This means that under the optimal 
stabilization conditions, interference suppresses not only 
ionization but also excitation of the level E2 , which can 
be seen experimentally also. 

In all pictures of this Section (Figs. 6-10) the calcu- 
lated values are plotted in their dependence on the in- 
tensity ratio X. To see such dependencies in experiments 
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one has to make, for example, a series of measurements at 
different values of the second-field intensity I2 at a given 
first-field intensity Ji . Another possible way of an exper- 
imental investigation is keeping the ratio x — I2/I1 con- 
stant and changing both intensities synchronously. Cal- 
culated for such a scheme of measurements, the proba- 
bility of ionization in its dependence on Ii oc I2 is given 
by the curves of Fig. 11. In this picture the intensity Ii 
is expressed in units of an arbitrary constant intensity Iq 
at a; = I2/I1 = 3. The detuning A and pulse duration r 
are taken to be equal to A = —200 x Iq and t — O.I/Iq 
(a) or r = I/Iq (b), where A, r, and Iq arc in atomic 
units. 

Normalization by an arbitrary constant Iq reflects the 
scaling effect described above in Section V. A possibility 
to choose any value of Iq indicates a large flexibility of the 
system under consideration with respect to a choice of of 
the light intensities and pulse durations. For example, if 
Iq = IQ-^a.u. (w 3 X IQi" W/cm^), the intensity I2 = 
3 X /i does not exceed 3 x 10^^ W/crn^ in the variation 
range of Ii at Fig. 11, which is low enough for no ATI 
effects to take place. And, under the best stabilization 
conditions, the detuning A and pulse duration are equal 
to A = -240 x Jo = -2.4 x 10""* a.u. - 0.065 eV <C uii^ 2 
and T — 1/Iq — lO'^a.u. ~ 3ps. 




Fig. 11. M,(7i/7o) at I2 = 3/1, S = -300 x (/0//1) and 
6» = 0.1 X (/1//0) (a) and h/Io (6). 

The curves of Fig. 11 are typical for the stabilization 
picture. With a growing light intensity, at first, the prob- 
ability of ionization grows (perturbation theory region), 
then falls, and this is the beginning of the stabilization 
window, and finally grows again, which corresponds to 
the break of stabilization. Stabilization and its break 
arise because, owing to the ac Stark shift and level mix- 
ing, with a growing light intensity the system comes to 
and, then, goes out of the resonance conditions, optimal 
for fuUy-dressed-level narrowing and interference stabi- 
lization. The curve b of Fig. 11 indicates an appearance 
of an additional region between the perturbation theory 
and stabilization zone at a sufficiently long pulse dura- 
tion T. This is the region of the total ionization of an 
atom, where Wi ^ I. This means that at some inter- 
mediate intensity light pulses provide a complete ioniza- 



tion and at a stronger field, owing to interference, ion- 
ization becomes rather small {wi ~ 0.2). In absolute val- 
ues, the minimal achievable probability in the stabiliza- 
tion region is somewhat lower in the case of short pulses 
(r — O.I/1//0) (the curve a of Fig. 11). But the degree 
of stabilization can be determined alternatively as the 
ratio of the maximal probability achievable in the region 
between perturbation-theory and stabilization regions to 
the minimal value of Wi in the stabilization window. In 
terms of such a definition, the degree of stabilization is 
much higher in the case of longer pulses [t — Ix/Iq) (the 
curve b of Fig. 11). Of course, a further increase of the 
pulse duration flattens the curve of Fig. 11 in the stabi- 
lization region and decreases the degree of stabilization. 
In this sense, the pulse duration chosen for the curve b of 
Fig. 11, T = Ii/Iq, is close to the optimal one. 

The picture of Fig. 12 characterizes spectral features 
of the residual probability to find an atom in its bound 
states at the best stabilization conditions of Fig. 11: 
I2/I1 = 3, h/Io = 1.2, 9 = 0.12 (a) and 1.2 (6) (which 
corresponds to r = O.I//0 and r = l/^o)- 
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Fig. 12. Wres{S) at I2 = 3/1, Ii = 1.2/0, = 0.12(a) and 
1.2 (fe). 

The curve a of this picture looks similar to the Fano 
curve of Fig. 1. This shows that at r = 0.12//o, still, 
the effect of stabilization under discussion can be inter- 
preted as a strong-field LICS. But in the case of longer 
pulses, T = I.2//0 (the curve b of Fig. 12) similarity with 
LICS practically disappears. The only reminder about a 
remote connection with the Fano curve is a slight asym- 
metry of the curve b of Fig. 12. Apart from this, the 
curve b describes the effect of interference stabilization 
in its pure form. 



C. Smooth envelope 

All the results described above were derived for pulses 
with a rectangular envelope. Usually, envelopes of short 
laser pulses are smooth. To consider such a more realis- 
tic situation, we have solved general equations (|ll|l with 
the sin^ pulse envelopes The results of such a so- 

lution are shown in Fig. 13, which is a direct analog of 
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Fig. 9. Again, a series of resonance-like curves describes 
the residual probability of finding an atom in its bound 
states Wres x) at various given values of the detuning 
S and the intensity ratio x considered as the indepen- 
dent variable. The residual probability maximized with 
respect to the detuning 6 is the envelope of the peaks of 
these curves. In Fig. 13 such a maximized probability 
is approximated by the functions Wres i^lpt (x) , x) , where 
^opti^) is empirically found linear function providing 
the best fitting to the peak envelope: 

6lft{x) =3.122 -73.3 X (49) 



w„{8, x), w„[8„^,(x), x] 




I I I 1_ 

0.5 1 1.5 2 



Fig. 13. A series of curves WresiSjX) at various given values 
of the detuning 5 and the function ■Wrea{So^jS{x),x) (a thick 
curve) for a He atom and the sin'^ pulse envelopes (Eq. 1381 . 

In Fig. 14 we plot the maximized residual probability 
to find He atoms in bound states calculated in the cases 
of rectangular (the curves 1 and 2) and sin^ (the curve 




0.3 - 



X 

0.2' ' ' ' ' ' — 

1 2 3 4 5 

Fig. 14. The functions <|f(^^;f(:E), a;) (1), w;ifiS^J;\{x),x) 
(2), and wlZ{Sopt{x) , x) (3) with Sl^fix) and J^j," (x), a; given 
by Eqs. I48II and 1491 . respectively. 

3) pulse envelopes at ^ = S^':f{x) (the curve 1) and 



6 = Sl^pf{x) 149() (the curves 2 and 3). Comparison of 
the curve 1 and 3 shows that a transition to a smooth 
envelope So{t) reduces a little bit the maximal achiev- 
able degree of stabilization compared to the rectangular- 
envelope case, but not too much (70-75 % instead of 
90%). Moreover, the curve 3 of Fig. 14 shows that in the 
case of a smooth pulse envelope the residual probability 
remains more or less stable in a rather large variation 
interval of the intensity ratio x, approximately from 0.5 
to 5 and more. This shows that the effect of stabilization 
is rather robust. 

Another interesting effect seen rather well from com- 
parison of the curves 2 and 3 of Fig. 14. These two curves 
are calculated at the coinciding dependencies of the de- 
tuning S on the intensity ratio parameter x, S = Sl^p^{x) 
(|49|l . As it's seen well from Fig. 14, at a; < 1 the curve 
3 goes above the curve 2. This means that at the same 
detunings the residual probability to find an atom in its 
bound states in the case of smooth envelope pulses ex- 
ceeds the same probability at a rectangular pulse enve- 
lope. In other words, in this range of the intensity ratio 
parameter x smoothing of the pulse envelope increases 
rather than reduces the degree of stabilization. This con- 
clusion follows directly from calculations though it looks 
counterintuitive and, in this sense, rather interesting. 

IX. HYDROGEN 

In a hydrogen atom, all the polarizability tensor com- 
ponents are known for the levels 2s and 5s and fre- 
quencies uji = 4.02 eV {XeCl laser) and W2 = 1-17 eV 
{Nd : Y AG laser) In atomic units they are given by 

ai(wi) = -45.56 + ^27.29, 01(^2) = 179.92, 
a2{uJi) = - 45.66 + i 1.78, a2(cj2) = - 513.76 + i 93.83 
and ai2 = 6.56 -I- i 50.60. 

(50) 

The data H50|) show, in particular, that Fj^^'/ri = 
(32 (wi)/q;"(u;i) ^ 6.5 X 10~^, which means that the as- 
sumption H18|) is pretty well satisfied. 

Rigorously, in a hydrogen, there are other levels (5d 
and 5/) of almost the same energy as 5s. Owing to the 
selection rules, the level 5/ is not connected either with 
2s or 5s levels by two-photon Raman-type transitions 
and, for this reason, can be ignored. As for the level 5d, 
in principle, it can participate in a scheme of two-photon 
Raman-type transitions under consideration. However, 
by ignoring at first this additional level let us consider a 
two-level 2s-5s scheme, analogous to that of the previous 
Section. _ 

For this system the coordinates of the (5 = point in 
the {x, 5} plane are given by 

xq = 0.272 and 5q = -47.2. (51) 

At (5 = (5o, the calculated relative width of the fully 
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dressed quasienergy levels g± H34() in their dependence 
on X are shown in Fig. 15. 

Compared with Fig. 5, the picture of Fig. 15 indicates 
the first well pronounced difference between Helium and 
Hydrogen. In the case of a Hydrogen the curves of widths 
of quasienergy levels vs. x have two avoided-crossing 
points whereas in the case of Helium such a situation 
never occurs and a,t S = Sq there are two real-crossing 
points (Fig. 5). 




0.1 0.2 0.3 0.4 0.5 0.6 



Fig. 15. The functions Pi (5o, x) 13311 

Another important difference concerns smooth pulse 
envelopes and the third level effect. To solve such a 
problem, we have to generalize Eqs. (|10|l and (|ll|l . In 
Eq. H1U|) . in accordance with the more general Eq. 
there appears an additional term C3(i) e*"^* ^3. Then, 
equations for Ci{t) {i = 1, 2, 3) take the form 

iCi - {Ei{t) + uji) Ci 
= -leio{t)e2o{t)(^ai2C2 + arse's), 

iC2-{E2{t)+uj2)C2 = -\eio{t)e2o{t)a2iCi 
ids - (Esit) + UJ2) = -i£io(t)e2oWai3Ci 

- 1 ("23 (C^l )£? (i) + «23 {UJ2)SI q (t)) C2 (t) 

where is given by the same Eq. (|12|) as Ei and E2 
with 03(0^1^2) and new off-diagonal elements of the po- 
larizability tensor given by '9'] 

asiuji) = -43.26 + i 0.42, 03(1^2) = -405.9 i 81.8, 

a23(wi) = 0.74 -I- i 0.21, a23(u'2) = 68.61 + i 21.63, 

and ai3 = 6.15 -f i 11.69. 

(53) 

Found from Eqs. (|52l) and (|35ll the residual probability 
to find an atom in its bound states is shown in Fig. 16 
for three different values of the pulse envelope smoothing 
factor a. 

Two rather interesting conclusions can be deduced 
from this picture. First, a strong smoothing of pulse 



envelopes decreases the peak value of the residual prob- 
ability to find an atom in bound states. At the chosen 
value of the detuning S = —530 in the case of pure sin^ 
pulses ['Wresix)]max IS almost twicc smaller then in the 
case of rectangular envelopes. The smoothing induced 
decrease of the residual probability in the case of Hydro- 
gen is much more significant than in the case of Helium. 




1 2 3 4 5 6 

Fig. 16. The function Wresix) in a three- level scheme at 
S = —530, 6 = 0.1, and the envelope smoothing parameter in 
ll^ a = 100, 10, and 0.1 (from top to bottom) . 

The second effect seen in the picture of Fig. 16 con- 
cerns the influence of the third level. This influence man- 
ifests itself in a shoulder on the curves Wres{x), but the 
third level is seen not to affect much the main maximum 
of the curves. 



X. CONCLUSION 

To summarize, we describe and discuss a scheme of in- 
teraction of atoms with radiation of two lasers. Intensity 
and pulse duration of lasers are assumed to be high and 
long enough to provide full ionization in the field of each 
of these two lasers alone if only atoms are prepared ini- 
tially at levels from which one-photon ionization can take 
place. We show that owing to interference effects under 
the conditions close to Raman-type resonance between 
some two selected atomic levels ionization of an atom 
experiencing a joint action of the field of two lasers can 
be significantly (up to 90 %) suppressed. Optimization 
of such a stabilization effect involves optimization with 
respect to the Raman-type resonance detuning and the 
ratio of the two laser intensities. Specific calculations are 
carried out for Hydrogen and Helium atoms for couples 
of atomic levels and laser frequencies at which informa- 
tion about the complex polarizability tensors involved is 
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available. Qualitatively, the results of calculations for 
Hydrogen and Helium appear to be very similar. This 
gives us a reason to think that the effect described is 
rather universal, and can occur also at other atoms, lev- 
els, and frequencies. The dependence of the effect on 
laser pulse shapes is investigated. It is shown that in the 
case of Helium atoms sensitivity of the results to a pulse 
shape is lower than in the case of Hydrogen atoms. In 
Helium, even in the case of smooth pulses, the degree of 
stabilization remains rather high (more than 70 %), and 
the effect exists at this level in a rather large range of 



the intensity ratio parameter x. The described scaling 
effect gives a possibility to select ranges of variation of 
the laser pulse peak intensities and pulse duration in a 
ranges most convenient for experimental observation. 
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